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Abstract. We derive new upper limits on the ellipticity 
of pulsars whose braking index has been measured, more 
tightening than those usually given, assuming that both 
' a gravitational torque and an electromagnetic one act on 
them. We show that the measured braking indexes of pul- 
sars are recovered in this model. We consider the electro- 
magnetic torque both constant and varying with time. At 
' the same time constraints on the pulsars initial period and 
on the amplitude of the gravitational waves emitted are 
obtained. 

Key words: magnetic fields - gravitational radiation - 
stars: pulsars: individual: (Crab, Vela, PSR 0540-69, PSR 
■ 1509-58) 
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1. Introduction 

From the timing measurements we deduce that radio pul- 
sars are subject to a systematic secular spin-down. The 
standard formula for pulsar spin-down is 



n = Kn n 



(i) 



The quantity being I the star momentum of inertia 
with respect to its rotation axis, is the torque acting on 
the st ar an d we will refer to K, following Allen & Hor- 
vath ( 1997 ), as the "torque function". K and the braking 
index n depend on the mechanism which is at work. For 
instance, n — 3 for pure magnetic dipole radiation and 
n = 5 for gravitational radiation. The braking index can 
be calculated, at least in principle, from the relation 



(2) 



which can be very easily derived from Eq. ([!]) . In practice 
only for four pulsars the measure of the braking index 
is not dominated by the timing noise (Lyne et al 
(Boyd et al. |1995| ),(Kaspi et al. |1994| ), (Lyne et al 



1993) 



1996) 



PSR 0540 - 69 : n = 2.28 ± 0.02 



PSR 1509 - 58 : n = 2.837 ± 0.001 



Vela : 



n = 1.4 ±0.2 



All these values are less than the canonical value, 3, which 
we expect for the emission of dipolar electromagnetic ra- 
diation. As we will see in the next section, different mech- 
anisms have been proposed to explain this discrepancy. 
On the other hand, when pulsars are considered as pos- 
sible sources of gravitational radiation, the upper limit 
on the amplitude of the waves emitted is calculated as- 
suming that all the spin-down is due to the emission of 
gravitational waves, e.g. (Hacnscl 1995] ), (Giazzotto et al. 



1997). This implies, see Eq.(14), a braking index equal to 



5. Clearly, this is not true, at least for the pulsars whose 
braking index has been measured. A much more realis- 
tic hypothesis is that in a pulsar both an electromagnetic 
and a gravitational torques work to produce the observed 
spin-down. We will show how from this assumption the 
observed braking indexes can be derived and, at the same 
time, how new more tightening upper limits on the pulsars 
ellipticity are obtained. 

The plan of the paper is the following. In Secj^ we give 
the main relations relative to the electromagnetic braking 
and shortly describe different proposed mechanisms which 
could produce the observed braking indexes. In Sec.^, af- 
ter the introduction of the relations describing the gravita- 
tional torque, we show how the observed braking indexes 
can be obtained combining the electromagnetic and gravi- 
tational torques. In Sec.£l| we calculate, in the framework of 
our model, the time evolution of various observable quan- 
tities, in particular the braking index n(t) and the angular 
velocity of the pulsar f2(£), and derive new upper limits on 
the ellipticity of pulsars and then, Sec.(||), on the ampli- 
tude of the gravitational signals emitted. Finally, in SecJ| 
conclusions are discussed. 
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2. The electromagnetic torque 



Assuming n em — const, from Eq.(M) we find 



As it is well known, the electromagnetic torque implies 
an evolution of the pulsar rotation frequency given by the 
relation 



(3) 

The electromagnetic index is n em = 3 for pure dipolar 
magnetic radiation (Ostriker & Gunn 1969| ), (Manchester 
& Taylor 1977), in which case 



2 



sin 



(a) 



(4) 



In this equation B is the strength of the magnetic field, 
R is the neutron star radius, a is the angle between the 
rotation and the magnetic axes, and I is the momentum 
of inertia with respect to the rotation axis. An electro- 
magnetic braking index different from 3 can be the conse- 
quence of pulsar winds, in which particles having angular 
momentum are accelerated away from the pulsar (Goldre- 
ich & Julian 1969), (Kaspi et al. 1994) or can be pro- 



duced if non dipolar components of the magnetic field are 
present, as a consequence o f stro ng magnetospheric cur- 
rents (Blandford & Romani 1988| ). Another possibility is 
that the magnetic moment, and then the torque applied 
to the star, varies in time (Blan dford & Romani 1988| ), 
(Cheng [1989|) . Muslimov & Page ( |1996| ) consider the time 
evolution of the surface magnetic field of the star, which, 
according to their model, is very low (10 8 -j- 10 9 G) for a 
newborn neutron star and then increases due to the ohmic 
diffusion of the initially trapped inner magnetic field. The 
resulting electromagnetic braking index is given by 



(5) 



where B sur f is the value of the surface magnetic field at 
the magnetic pole. Clearly, if B sur f > 0, n ern < 3. Other 
models based on the secular variation of the magnetic 
fiel d have been developed , for instance, by Blandford et 
al. ( |1983|) and by Camilo (|l996|) . A different kind of mod- 
els has been pro posed by Allen & Horvath ( 1997 ) and by 
Link & Epstein ( 1997 ), based on the growth of the angle 
a between the magnetic moment and the rotation axis of 
the star. In this case the electromagnetic braking index is 




= 3 



tana \ft 



n 



(6) 



In all these cases the general expression for the observed 
electromagnetic braking index is simply 



n eni — <jH — X 



n 



k, 



(7) 



which is easily obtained deriving Eq.(g) with respect to 
time. From Eq.fl) we find 



K„ 



1 



B sin a 



2£?sin a + B sin 2a a — B sin a- 



(8) 



.,0 x 



(9) 



where K em and fl are the electromagnetic torque func- 
tion and the angular velocity at time t — 0, i.e. now. Then, 
the rate of variation of the angular velocity can be written 

as 

ftem = K em ^ e ff£l n ' m (10) 



where K em<e ff — n n^™'° 3 is a constant, which is formally 

equivalent to Eq.(||). Integrating Eq.(||) we obtain the pul- 
sar "characteristic" time, due to the electromagnetic brak- 
ing which is: 



n 



(1 - n em )fi 



(11) 



being f2i the initial angular velocity of the pulsar. Eq.([LlJ) 
would give the true age if only the electromagnetic emis- 
sion was responsible for the pulsar spin-down. 

From Eqs.(|]||), tne strength of the magnetic field, for 



3, is 



B sin a ~ 3.2 • 10 19 \/ PP em G 



(12) 



where a is the angle between the magnetic and the rota- 
tion axis, P em is the rate of variation of the pulsar period 
due to the electromagnetic torque, and we have assumed 
I = 10 38 kg m 2 . 



3. Electromagnetic plus gravitational torque 

In this section we show that the observed braking indexes 
can result from the combined action, in a pulsar, of the 
magnetic and gravitational torques. We assume that the 
emission of electromagnetic and gravitational radiation 
are the only mechanisms acting in the pulsar. The ob- 
served spin-down rate can be expressed as 



n = n e 



n 



where tl em is given by Eq.(||), or Eq. (|To|) , while 
32 G 



n 



5 c 5 



-ie 2 n 5 = K n ,„n 5 



(13) 



(14) 



where e is the ellipticity of the pulsar. We note that, in 
analogy with Eq.jTl|), the pulsar "characteristic" age, due 
to the gravitational braking, is 



go 
4£7 



1 



ft. 



(15) 



In the following, we will consider the torque functions K e 



and K gw constant in time. 
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Fig. 1. Y = as a function of the braking index n for 
different values of the electromagnetic braking index n em . 
n em varies in the range [1, 3] with step 0.2. 



The braking index n can be obtained differentiating 
Eq.(|l3|), and using Eqs.(|Q: 



ft x ft n em + 5Y 

n = : = 

f]2 1+Y 
where we have introduced 



gw 



K„ 



K, 



(16) 



(17) 



This quantity is the ratio between the gravitational and 
the electromagnetic contributions to the pulsar spin-down 
rate and we will see in the next section that its present 
values are always less than f . Note that, in terms of f2 and 
its derivative, 



Inverting Eq.(16) 



Y(n) 



(18) 



(19) 



which implies n em < n < 5. In Fig. (|T|) we have plotted 
the function Y(n) for different values of n em in the range 
[1, 3]. As we expect, for each n em , Y(n) = when n = n em 
(i.e. no gravitational torque) while Y(n) — > +oo for n — > 
5 (i.e. negligible electromagnetic torque). We see that a 
braking index less than 3 can be obtained combining the 
effects of the electromagnetic and gravitational torques. 
For instance, for the Crab pulsar we have n = 2.51 and 
if we assume, say, n em — 2.4, then it immediately follows 
Y ~ 0.04, that is the gravitational radiation contribution 
to the spin-down rate is about 4% of the electromagnetic 
one. On the other hand, if n ern = 1, then Y ~ 0.6. 



The spin-down of a pulsar in which both the electro- 
magnetic and gravitational torques are acting, Eq.(p"3"|), 
can be expressed as 



n = K eff (Q) x Q n{n) 



(20) 



where both the "effective" K e f f and n are function of the 
angular velocity f2. As a consequence of its dependence on 
f2, the braking index n is not constant in time. In the next 
section we will derive the time evolution of n. In all subse- 
quent calculations we will assume the pulsar momentum of 
inertia being equal to the canonical value / = I0 38 kg m 2 . 
From Eq. (fL4J) we can express the ellipticity of the pul- 



e = 1.9 x I0 5 \/P 3 P ( 



(21) 

Combining Eqs. ( |14| , |17|) we can re- write Eq.(^l|) as a func- 
tion of the observed pulsar period and its derivative, and 
of the ratio Y between the gravitational and electromag- 
netic spin-down rates: 



e = I.9xl0 5 WP 3 P- 



Y(n) 



= 7.55x10 




Y(n) 



1 + Y(n) " y n 5 1 + Y(n 

Eq. (^||) plays a basic role in the determination of the up 
per limits of the pulsars ellipticity. 

4. Pulsar parameters evolution 

Differentiating Eq. (|l7j) we easily find 



(22) 



Y = (5- n em )Y^ 
and then 



where Yq = 



(23) 



(24) 



5-no 



being hq the observed braking in- 
dex (i.e. its present value). We use this solution to solve 
the equation for the time evolution of the pulsar angular 
velocity. To this purpose, let us write Eq.([13|) in the form 



-l^L gqjj 



Y(Q) J 



(25) 



We numerically solve this differential equation, fixing the 
pulsar age t age and finding a solution depending on three 
unknowns: the ellipticity e (contained in K gw ), n em (on 
which Y(Q) depends) and the initial angular velocity Qi 
(Q i = Q(t age )). These quantities must be chosen in order 
to verify two conditions: those given by f2(i = 0) = Qo and 
by Eq.(p3), calculated for the values of the parameters at 
t = 0Q: 



7.55 x 10 6 i 



|n | Y 



1 + Y 



(26) 



In fact, as we assume K gw = const, the value of e coming 
from Eq.(^) is not dependent on the particular time t consid- 
ered. 
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This gives us a range of acceptable values for the un- 
knowns. The solution we are searching for is that corre- 
sponding to the smallest possible value of electromagnetic 
braking index n em , that is to the greatest value of the el- 
lipticity eQ. Once we have determined the function fi(t) 
for a given pulsar, we can use Eq.([l6|) and find the time 
evolution of the braking index n(t). 

In some cases, we can further reduce the range of vari- 
ation of n em considering the pulsar age, in the following 
way. Let us consider the ratio R between fl em and the 
total spin-down rate, given by Eq.(^), or equivalently, by 



n = K em n n ™ (i + y(si)) 

We have 



i 



i + y(n) 



(27) 



(28) 



5=no_ < j_ Thc 

5 — Hem — 



This quantity varies monotonically between the two limits 

Rinf = R{&i) anQl Rsup — R{Qo) 
condition Vt ern < R sup Q implies 

^aae 



Rsup ^* 



(29) 



with r em given by Eq.(pd|) . We choose the range of n em in 



order to satisfy Eq.(|29|). In this way, with the exception 
of the Vela pulsar, we can increase the lower limit of n em 
then reducing the upper value of e. 

In the following we consider individually the four pul- 
sars with measured braking index, applying to each of 
them the procedure described. 

4.1. Crab 

The Crab pulsar present parameters are n n = 2.51 ± 0.01, 



n = 

is t a 



188.12 s~ x and fin 



-2.42 x 10- y s -2 , while its age 



t 945 yr. The initial rotation period of the Crab is 
not exactly known but a value R t £ [10, 20] ms is usually 
assumed. By solving Eq.(p5[), under the various conditions 
previously described, we find n em > 1.7, e < 3.8 x 10~ 4 , 
Cli < 394 s _1 (Pi > 16 ms). A more tightening limit can 
be obtained considering that the observed braking index 
appears to be constant within 0.5% in more than twenty 
years of obs ervati ons, thus implying |n| < 5 x 10~ 4 yr^ 1 
(Lyne et al. 1993 ). This relation is satisfied if n em > 2.05 
and correspondingly we find e < 3 • 10~ 4 , 2.5 times lower 
than thc maximum ellipticity calculated under the hy- 
pothesis that the pulsar spin-down is due only to the emis- 
sio n of g ravitational waves (Haensel 1995 ), (Giazzotto et 
al. |1997|) , Q; < 364 s^ 1 (P % > 17.2 ms) and Y < 0.18. 



age is not well determined but estimates, based on the 
position offset and the proper motion of the pulsar, give 
values in the range [18,31] x 10 3 yr (Aschenbach et al. 



|1995 ). Assuming t age = 1.8 x 10 4 yr, we have n em > 1.27, 



e < 3.3 x 10 , about 5.5 times lower than the limit relative 
to the case of gravitational torque alone, fii < 391 s~ 4 
(Pi > 16.1 ms) and Y < 0.036. If t age = 3.1 x 10 4 yr 
we find much more tightening limits: n em > 1.399, e < 
3 x 10~ 5 , about 60 times lower than the limit relative 
to the case of gravitational braking alone, together with 
Cli < 594 s- 1 (Pi > 10.6 ms) and Y Q < 2.78 • 10~ 4 . 



4.3. PSR 0540-69 

The pulsar PSR 0540 - 69 is characterized by n = 2.28 ± 
0.02, fi = 124.66 s" 1 and O = -1.185 x 10~ 9 s" 2 . We 
do not know exactly its age which, however, should be of 
about 900 yr (Muslimov & Page |1996[ ). We find n em > 
2.21, e < 2.4 x 10~ 4 , which is 6.3 times lower than the 
"gravitational" limit, fii < 175 s _1 (Pi > 36 ms) and 
Y < 0.026. 



.4. PSR 1509-58 



2.837 ± 



10 s- 2 . If the 



The pulsar PSR 1509 — 58 parameters are: no 
0.001, n = 41.7 s- 1 and 57 = -4.25 x 10 
pulsar birth is identified with the supernova event AD 185, 
then its age is t age ~ 1800 yr (Thorsett 1992), (Muslimov 
& Page 1996). On the other hand its "electromagnetic" 



age is T em < 1700 yr. The origin of this discrepancy has 
not been understood yet. So, our model cannot be applied 
because the pulsar age resulting from the action of both 
the gravitational and electromagnetic torques is always 
less than r em and cannot match t age . 

We have plotted the time evolution of thc pulsars an- 
gular veocity f2(f) in Fig.(|2|) and that of the braking in- 
dex n(t) in Fig.(||). In all cases n(t) — * n em asymptoti- 
cally, because the gravitational contribution to the spin- 
down rate decreases faster than the electromagnetic one, 
see Eq.(|r7j).We note that for the Vela pulsar, if t age = 
1.8 • 10 4 yr, the initial value of the braking index n could 
be as high as 4.8. This implies that Yi = V(f^) < 22, so 
that at the beginning the spin-down could have been dom- 
inated by the gravitational braking. From the figure we see 
that Y would have become less than 1 about 3000 yr af- 
ter the birth. Also for the Crab the gravitational braking 
could have been dominant at the birth (Yi < 1.28). 

Our results, together with some pulsars parameters, 
are summarized in Tab.rtl). 



4.2. Vela 

The Vela pulsar is characterized, at present, by no — 1.4± 
0.2, n a = 70.36 s- 1 and Cl = -9.85 x lO" 11 s~ 2 . Its 

2 From Eq.(^), the greater is the difference between n em 
and no and the greater is the ellipticity. 



5. Emission of gravitational waves 

From the calculated upper limits on the ellipticity of pul- 
sars we can very easily find new upper limits also on the 
amplitude of the gravitational signals emitted by them. 
For a deformed neutron star rotating around a principal 
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Fig. 2. The angular velocity fl(t) is plotted as a function 
of time, assuming n em equal to its lower limit. For the 
Vela pulsar age both the lower and the upper limit have 
been considered. 
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Fig. 3. The braking index n is plotted as a function of 
time, assuming n em equal to its lower limit. For the Vela 
pulsar age both the lower and the upper limit have been 
considered. 

axis of inertia, the waveforms corresponding to the two in- 
dependent polarization states are, in the quadrupole for- 
malism, 



h+ = yCos2m(l + cos 2 
h x = ho sin 2Vtt cos 6 



(30) 
(31) 



where 9 is the angle between the rotation axis of the star 
and the line of sight and the amplitude ho is given by 

I \ /10 kyc 



h Q = 1.05 x 10" 



10 3S kg m 2 



io- 



r)( 



100 Hz 



(32) 



being v = — the frequency of the wave, equal to two times 
the rotation frequency of the star. For completeness, we 
remember that the power emitted through gravitational 
waves is 



(33) 



dE = 32 G /W 
dt 5 c 5 

In Tab. flj]) we have reported, among the other quantities, 
the maximum amplitude of the gravitational waves emit- 
ted from the pulsars with known braking index, calculated 
by the use of Eq. (|3^) , assuming I — 10 38 kg m 2 . 

6. Conclusions 

Usually, the upper limits on pulsars ellipticity are calcu- 
lated assuming that pulsar spin-down is due only to the 
emission of gravitational waves. It is well known that this 
assumption contradicts the measure of the braking index, 
which has been obtained for a small sample of pulsars, 
with values always less than 3. 

In this paper we show that considering the simulta- 
neous action, in a pulsar, of a gravitational torque and 
an electromagnetic one, the observed braking indexes can 
be recovered. The gravitational torque is always a small 
correction to the total torque acting on the pulsar. A con- 
sequence of the model is that the braking index is not con- 
stant in time. We calculate the time evolution of the pul- 
sars angular velocities and braking indexes and find more 
tightening upper limits to their ellipticity. These new lim- 
its are between 2.5 and 60 times lower than the old one. 
New upper limits on the amplitude of the gravitational 
waves emitted immediately follow. 

Acknowledgements. I would like to thank the anonymous ref- 
eree for his remarks and his helpful suggestions. 
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r [kpc] 


Po [ms] 


v [Hz] 


Pi,min [?7is] 


Y(),rnax 


£max 


hmax 


Crab 


2 


33.4 


59.9 


17.2 


0.18 


3 x 10~ 4 


5.5 • 10~ 2!) 


Vela 


0.5 


89.3 


22.4 


i6.r 


0.036 


3.3 x 10~ 4 


3.5 x IQ~ Z0 










10.6" 


0.00028 


3 x 10~ b 


3.2 x 10~ 2b 


PSR 0540 - 69 


49.4 


50.4 


39.7 


36 


0.026 


2.4 x 10" 4 


8.1 x lO" 27 
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emitted. For the Vela pulsar two ages, t age = 1.8 x 10 4 yr (a) and t age = 3.1 x 10 4 yr (b) have been considered. 
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